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We consider the non-equilibrium steady states of a driven charge density wave in the presence of 
impurities and noise. In three dimensions at strong drive, a true dynamical phase transition into a 
temporally periodic state with quasi-long-range translational order is predicted. In two dimensions, 
impurity induced phase slips are argued to destroy the periodic "moving solid" phase. Implications 
for narrow band noise measurements and relevance to other driven periodic media, e.g. vortex 
lattices, are discussed. 
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PACS numbers: 71.45.Lr,72. 70. +m, 74.60. Ge 

The influence of quenched impurities on a periodic 
medium can lead to very rich physics. Examples include 
charge density wave (CDW) systems and the mixed 
state of type II superconductors § , in which the vortices 
form a periodic lattice. In both these cases, it has been 
argued that the impurities ultimately destroy the long- 
ranged periodicity and pin the periodic medium. How- 
ever, with an applied force, provided by an electric field 
or current, the periodic structure de-pins and becomes 
mobile. Once in motion, the impurities are less effective 
at destroying the periodicity. Indeed, recent experiments 
on YBaCuO have shown evidence for a first order melt- 
ing transition of the moving vortex lattice ||. Vinokur 
and Koshelev have interpreted this experiment in terms 
of a true non-equilibrium phase transition, from a moving 
liquid phase to a "moving solid phase" Q . 

This experiment raises a number of questions about 
the non-equilibrium steady states of such noisy driven 
systems with impurities. The most basic concerns the 
very existence of a "moving solid" phase. A solid in equi- 
librium is usually characterized by the presence of long- 
ranged crystalline correlations (Bragg peaks). But other 
criteria also suffice, such as the presence of a non-zero 
shear modulus, or the absence of unbound dislocation 
loops. Under what circumstances, if any, is it possible to 
have a true moving solid, separate from a driven liquid 
with plastic flow? If the moving solid phase is possible, 
what are it's characteristics and experimental signatures? 

In this letter we attempt to answer these questions, 
focusing for simplicity on the CDW. Many of our con- 
clusions, however, apply also to the driven vortex lattice. 
In two dimensions, we find that a moving solid phase 
driven through impurities is always unstable to a prolif- 
eration of dislocations. The system becomes equivalent 
(in symmetry) to a driven liquid. In three dimensions, 
a moving solid phase appears to be stable at large ve- 
locities, as illustrated in the schematic phase diagram, 
Fig.l. However, the solid phase does not have true long- 
ranged positional correlations (LRO), as in an equilib- 
rium (3d) crystal. Rather, algebraic power law positional 
correlations are predicted, as in a 2d equilibrium crys- 
tal. Likewise, unbound dislocation loops are absent in 
the moving solid. Despite the absence of spatial LRO, 
the moving solid phase is periodic in time - and hence 



has long-ranged temporal correlations. The experimental 
signature of such a periodic state is narrow band noise, 
at the "washboard" frequency ^Bl. 
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FIG. 1. Schematic 
phase diagram for the three-dimensional CDW. Hashed lines 
indicate (possibly sharp) crossovers. 

Upon inclusion of thermal effects or phase slips, the 
CDW de-pinning transition |(| is predicted to be rounded 
0, becoming a crossover (see Fig.l). For electric fields, 
E, above this crossover, the CDW is in a plastic flow 
regime. With increasing E, we predict a true phase tran- 
sition into a temporally periodic "moving solid" phase. 
For the CDW, this dynamical transition is likely to be 
continuous. As shown below, scaling arguments then pre- 
dict NBN characteristics near the transition. The phase 
diagram for the driven vortex lattice should be similar, 
with current replacing electric field. However, in this case 
the transition into the moving solid phase is likely to be 
first order, at least in the large current limit. 

Charge density waves tend to form in very anisotropic 
metals, consisting of weakly coupled metallic chains [||. 
The electronic density in a CDW has a periodic modula- 
tion along the chain (x-)dircction: 



p(x) = po + Retpe 



(1) 



with kp the in-chain Fermi wave vector. Long-ranged 
order of the CDW is manifest in the complex order pa- 



rameter field, "0(x) = pie 
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In the absence of impurities and an applied electric 
field, the CDW exhibits long-ranged order in the pair 
correlation function: 



G(x) =< V*(x,t)V(0,t) > t , 



(2) 



with G(x — ► oo) 7^ 0. Here the subscript t denotes a 
time average (equivalently an ensemble average in equi- 
librium). Lee and Rice have argued that quenched im- 
purities destroy the LRO of G, for physical dimensions, 
d < 4 ||. However, when the CDW is driven and mov- 
ing, the Lee-Rice argument is not valid, and LRO of G 
is not precluded (but see below). 

In the moving non-equilibrium steady state, temporal 
correlations in ip also serve to characterize the CDW or- 
der. Consider the pair correlation function: 



Git) =< V*(x ) t)V(x,0) > a 



(3) 



where the subscript x denotes a spatial average. Tem- 
poral LRO is signaled by a periodic, and non-decaying, 
behavior: C(t + to) = C(t), with to the "washboard' 
period. Middleton has shown that for a large class of 
dynamical models, which exclude phase slip and ther- 
mal noise, the steady state is (microscopically) periodic 
j|. With noise microscopic periodicity is destroyed, but 
the statistical correlation function C(t) is still periodic. 
However, when phase slip is allowed, the robustness of 
the periodic state is much less clear, as we discuss below. 

To proceed, we assume local CDW order, and con- 
struct a long wavelength description in terms of ijj. In 
the absence of phase slips, as we consider first, ampli- 
tude fluctuations can be ignored, and the dynamics in- 
volves the phase field, (f>. A common starting point is the 
Fukuyama-Lee-Rice (FLR) model [pi, with equation of 
motion: 



d t cj) = DV 2 (f> + V(x) sin(2fc F a; + (f) + v2k F . 



(4) 



The spatial coordinates transverse to the chains have 
been re-scaled to give an isotropic diffusion constant 
D = TVp, with scattering time r and Fermi velocity v F . 
The second term on the right side represents the effect of 
quenched random impurities. The last term is present in 
an applied electric field E, with velocity v = (eT/m)E. 
This term can be shifted away, 4> ~ > & + v2kpt, reduc- 
ing the FLR equation to an equilibrium form: dt<p = 
—SH/6(f), with Hamiltonian H. However, there are ad- 
ditional terms that can, and should, be added to FLR, 
which arc manifestly non-equilibrium, as we now discuss. 

The most important such term is d x (f>, allowed by sym- 
metry once the CDW is in motion along the x-direction. 
The magnitude of this term can be estimated by bal- 
ancing forces on a single wavelength, A, of one chain. 
In an electric field E, the force is 2eE, accelerating (lo- 
cally) the phase field: 2md 2 (A0 / '2ir) = 2eE, or equiv- 
alently, rdfcj) = v2tt/X. In the presence of distortions 
in <p(x), the (local) wavelength of the CDW is modified: 
2ir/\ — * 2fcp + d x 4>. Thus the last term in FLR should be 
replaced by: 2fcp — * 2/cp + d x <fi. (Note that the inertial 
term has been dropped.) 



In general there are other missing terms, for example 
of the KPZ form, (d x <j)) 2 (il]] . This term, involving more 
gradients and powers of (f> is less relevant than d x <j>. In 
the following we drop this term, although it can play an 
important role jl^ ]. 

In the moving state, it is tempting to argue that the 
impurity term, V(x), will average to zero at long times. 
However, as pointed out by Coppersmith 0, it is not 
legitimate to ignore completely the effect of impurities. 
In particular, the impurities modify the local mobility, 
/i, of the phase field. (This can be seen explicitly via a 
high velocity expansion.) The random V term can then 
be replaced by 8fi(x)E2k-p = F(x), where S/i denotes the 
fluctuating part of the mobility. We take F Gaussian 
with F = and F(x)F(0) = g5(x). 

We thereby arrive at a generalization of FLR: 

d t( j) = D\7 2 <t) + v(2k F + d x <p) + F(x) + »7(x, t). (5) 

The stochastic noise term is assumed to be Gaussian with 
(r)) = and (r?(x, t)r;(x', t')) = k B T5(x - x')S(t - t'). 

Finally, we modify the model to allow for phase slip 
processes. One way to incorporate phase slip is to put the 
model on a lattice and replace d x (j) — > a~ 1 sin((f>(x + a) — 
4>(x)), etc. Alternatively, amplitude fluctuations can be 
included, using field ip. The appropriate soft-spin model, 
which reduces to Eq.|^ in the spin-wave limit, is, 



d t ip = [-DV 2 + vd x + M + r(x) + ioj Q + iF(x)]tp 
- U ^| 2 +£(x,i), 



(6) 



with the definition uio — v2kp. Here M is a "mass" term, 
which controls the magnitude of the order parameter, 
and £(x, t) is a complex stochastic noise term. We have 
also included a spatially random component to the mass, 
denoted r(x), which we take to be a zero-mean Gaussian 
random variable with r(x)r(0) = A<5(x). 

First consider the system with phase slips suppressed. 
In this spin- wave limit, EqJ^ is linear in </> and can be 
solved via Fourier transforms: 



F(P) 



Dp 2 + ivp : 



-2it5(uj) 



iuj + Dp 2 + ivp x ' 



(J) 



where <fi(t) — <fi(t) — u>ot. The first term in Eq.0 repre- 
sents a static distortion in 4>{x) induced by the random 
mobility, while the second gives "noisy" dynamical fluc- 
tuations around this mean. The static phase variations 
diverge algebraically with system size L for d < 3, lead- 
ing to (stretched) exponential decay of G(x). Thus, even 
without phase slips, a 2d driven CDW lacks translational 
LRO. For the 3d case, EqJ^ gives G sw (x) ~ x~ v , corre- 
sponding to power law peaks in the static structure func- 
tion and translational quasi-LRO (QLRO). Notice that 
the presence of the non-zero vd x (f> term in Eq.|| is critical 
here. Indeed, for a driven periodic system with reflection 
symmetry, v = 0, Eq.|^ implies exponential decay of G(x) 
for all d < 4. Because the disorder term in Eq.|?] is static, 
the dynamical properties are determined by the thermal 



2 



noise term, with C sw {uj) ~ 5(ui — ojq) for all d > 2, in- 
dicating temporal LRO of C(t). In 2d spin- waves imply 
temporal QLRO for C(t). 

We now address the stability of these spin- wave results 
in the presence of phase slips. In equilibrium, v = F = 0, 
Eq.|| describes an XY model with relaxational dynamics. 
In this case the unbinding of topological defects (i.e.. vor- 
tices) coincides with the loss of translational LRO due to 
spin wave fluctuations. For d > 2, the vortices form d— 2 
dimensional subspaces (lines in 3d). With a core energy 
growing with size L as L d ~ 2 , they are bound at low tem- 
peratures. In equilibrium, 2d is marginal both for spin 
waves, which give QLRO, and for vortex unbinding. But 
in the non-equilibrium case of interest, the unbinding of 
phase slips and vortices needs to be re-addressed. 

For simplicity, consider first the case v — 0. Then 
Eq.[|can actually be cast into an equilibrium form, d t <j) = 
—SE c s/S<j) + T], with the proviso that the "energy" E e s = 
f{D\ V0| 2 /2 — F(f>} is a multi-valued (i.e. non-periodic) 
function of the phase. The situation is analogous to the 
"tilted washboard" model of Josephson junctions, except 
that the tilt is here a random function of position. 

It is clear that spin-wave conformations of the phase 
are highly constrained. Imagine sub-dividing the system 
into regions of linear size L. Each such region experi- 
ences a net random torque / d d rF(r) of order ±y^gL d . 
The torque in neighboring regions is generally different, 
so that the local phases are pushed at different rates. In 
the absence of phase slips, however, all regions must ro- 
tate synchronously or build up enormous strains. Eq.|?] 
describes the resulting steady state in which the strains 
increase to counteract the non-uniform applied torques. 

Once phase slips are allowed, however, such a situation 
is clearly metastable. If the net torque in a particular re- 
gion is positive, then the energy of the spin-wave state is 
lowered simply by increasing all the phases in the region 
L by 2tt. This decreases the random energy but does not 
alter the strain energy (which is now periodic in 4>). For 
finite L and non-zero T, this process will therefore occur 
with an activated rate 1/r ~ exp(— U/k&T), where U is 
the energy barrier for the phase slip process. 

The energy U is estimated from the elastic energy mid- 
way through the process, i.e. when there exist phase 
shifts of order ir on scale L. Adding the elastic and ran- 
dom contributions to the energy gives U(L) ~ DL d ~ 2 — 
\J gL d . A more microscopic picture is that of vortex nu- 
cleation. The phase slip is achieved by nucleating a small 
neutral topological defect (vortex-antivortex pair or vor- 
tex loop in d = 2,3), which expands and slips over the re- 
gion, annihilating again on the opposite side. The elastic 
contribution to the barrier energy is just the binding po- 
tential of the defect, Vdefect(-k) ~ DL d ~ 2 , up to possible 
lnL dependence. For d > 4 and small g, U(L) is positive. 
Moreover, since U(L) grows with large L large phase slips 
are exponentially suppressed, indicating stability of the 
spin-wave phase for d > 4. For d < 4, however, U(L) 
becomes negative for L > L c ~ (D 2 / g) 1 ^ 4i ~ d \ On scales 
much bigger than L c it is then inconsistent to assume a 



well defined average phase, because of the proliferation of 
vortex pairs/loops on smaller scales. The relaxation time 
for the phase slips on scale L c is r c ~ exp[Z?LJ?~ 2 /ksT]. 
Beyond this time scale, phases separated by distance 
large compared to L c will become de-phased, destroying 
the temporal LRO of the spin-wave state. 

For i; / 0, the argument is trickier. First, transform 
to the moving frame via x — > x — vt, which removes 
the vd x <j) term in Eqj^. The elastic force is invariant 
under such a transformation, but F (r) — > F (x — vt, rj_), 
so that the random torque field appears to move with 
velocity v. Again dividing the system into regions of size 
L, we see that a statistically uncorrelated realization of 
F moves into a particular region in a time to = L/v. 
For large L this decorrelation time to is much smaller 
than the typical diffusive time for phase changes, t$ ~ 
L 2 /D. Only on time scales longer than t<p can a phase 
change take advantage of the random torques spread out 
over the entire region. The random energy gained is thus 
averaged over t^/to realizations of the torques, leading to 
a net torque on the entire region of F net ~ \fgfvL ( d_1 )/ 2 . 
In this case the characteristic energy balance is U(L) ~ 
DL d ~ 2 — \JgJvL^ d ~ x i I" 2, , and phase slips proliferate for 
d > 3. Directly in d = 3, these naive arguments suggest 
a transition between an ordered state for v > v c and a 
disordered state for v < v c , with v c ~ g/D 2 . 

The above arguments are consistent with the spin- 
wave calculations. As in equilibrium, they suggest that 
vortex unbinding coincides with the loss of translational 
LRO due to spin-wave variations. Further support for 
this conclusion follows from an analysis of the soft spin 
model, Eq|j, as we now describe. 

In the absence of randomness, the soft spin model con- 
tains two phases for d > 2. Fluctuation effects near the 
transition, negligible above d = 4, can be studied for 
small e = 4 — d via the renormalization group (RG) . The 
RG including r(x) has been studied in the context of 



random-bond XY magnets |13|; we generalize this calcu- 
lation to include F in the dynamics (for v = 0). After 
transforming ip — > e ltJ °*'i/', we employ standard dynamical 
RG methods. The resulting differential RG flow equa- 
tions to quadratic order are, 

u = u(e — 5m + 6A), 
A = A(e - 4m + 4A - 2g) + 2g 2 , 
g=g(e-2g + 6A), (8) 

with u — du/d\nb, etc., where b — e e is the rescal- 
ing factor. A simple analysis shows that the Gaussian 
(u = A = g = 0), pure (A = g = 0) and dirty equilib- 
rium (g = 0) fixed points are all unstable, and further, 
that no other fixed points exist. Instead the couplings 
diverge as £ — > oo. In particular, u — > +oo (not — oo), so 
the instability does not appear to indicate a fluctuation 
induced first order transition. Instead, the strong diver- 
gence of the disorder strengths g and A are consistent 
with the scenario that the ordered phase is absent. 

For v 7^ 0, changing to co-moving coordinates x — > 
x — vt reduces Eq.|| to the previous case but with r(x) — > 
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r(x — vt,x±) and F(x) — > F(x — vt,x±). Because 
z = 2 + 0(e) at the pure XY fixed point, the weaker 
x dependence of r and F may be ignored. Perturbations 
of the form r(— vt,x±)tp and iF(—vt,x±)il> are strongly 
irrelevant near d = 4, consistent with the spin- wave anal- 
ysis and our earlier scaling arguments which gave d=3 
as the lower critical dimension for the ordered phase. To 
study d=3 in the soft spin representation, we employ non- 
perturbative techniques. We therefore consider a gen- 
eralized model containing N complex fields ipi obeying 
Eq.||, but with \tp\ 2 — > J2i IV^| 2 - We analyze the stability 
of the pure N — oo fixed point to the random pertur- 
bations. From scaling the singular part of the mean en- 
ergy density varies as (ip 2 ) = £, 1 ^^ d f(A^ VA ,gS, yg ), where 
£ ~ M~ x l v is the (pure) correlation length, and y& and 
y g are the RG eigenvalues of A and g, respectively. At 
N = oo, 1/v — d = —2. Differentiation implies A<\ = 

dAWjU, g =o ~ e yA ~ 2 , A g = d g W)W g =» r ?°- 2 . 

These quantities are computed at N = oo using saddle 
point techniques and the Martin-Siggia-Rose dynamical 
formalism |l4| . We find j/a = d — 5 and y g = 3 — d. (As 
a check we also considered the case v = 0, and found 
yA = d — 4 and y g — 4 — d, in agreement with the usual 
Harris criterion and the e-expansion, Eq.||.) Thus for 
d < 3, the equilibrium critical point is unstable to ran- 
dom F, consistent with the absence of an ordered phase. 

Our predictions for the 3d phase diagram are summa- 
rized in Fig.l. Upon lowering the temperature at weak 
drive, E, substantial CDW amplitude develops at the 
mean- field transition temperature T c q. At long distances 
and times, however, both G(x) and C(t) decay exponen- 
tially to zero. With increasing drive, the CDW undergoes 
a sharp non-equilibrium phase transition into an ordered 
"periodic state" with spatial QLRO and temporal LRO. 
Our arguments strongly suggest that for 2d CDW sys- 
tems, the ordered phase is absent. 

Experimentally, temporal LRO in the solid phase 
manifests itself in NBN. Consider current fluctuations 
in the presence of a fixed bias voltage; other set-ups 
are qualitatively similar. The current density j(x) = 
(enj_/7r)9t^[l + (/9i/2fcp)cos(2fcpa: + 0)], where n± is the 
areal chain density. In a sample of cross-sectional area A, 
the instantaneous CDW current through the plane x = 
is I x (t) = J A d 2 x±j(x = 0, Xj_, t). The oscillatory part of 
the NBN correlator, S(t) = (J* (i)/ x (0)) is thus 

S a (t)^I 2 [ Re{ e ^ t (V(xi,i)V'*(x' ± ,0))}, (9) 

where Iq = (enj_wo) / '(2v / 2fcF'?r) • We consider this quan- 
tity in the bulk, and expect that measured current fluc- 
tuations (in the external leads) exhibit proportional be- 
havior. Temporal LRO in the solid phase therefore im- 
plies a sharp (resolution limited) delta function peak in 
S(w) ~ A 2 ~ ti / 2 5(lj — uj ). Deep in the liquid phase, ip 
correlations are short-range in space and time, which 
gives the mean- field result S(u) ~ AVt/(VL 2 + uj ) 2 ). 
Near the transition field E C (T), provided the transition 
is continuous, we expect a scaling form S±(u>,SE) ~ 



\5E\ a f±[(uj - uj )/\5E\ zv ,A\5E\ 2 % where z and v are 
the dynamical and correlation length exponents, a is an 
additional scaling exponent, and SE = E — E c . Match- 
ing to the infinite area limit implies that the amplitude 
of the delta function frequency peak for E > E c van- 
ishes as \SE\ a+ '- 4:+z ^ n)l/ . For E < E c , matching im- 
plies the (generally non-Lorentzian) line-shape S(u>) ~ 
A\5E\ a+2w s[(uj - w Q )/\8E\ zv \. In two dimensions, S{u) 
has an intrinsic width for all fields and temperatures. 

Although the discussion has focused on CDWs, most 
of the ideas employed here apply to more general peri- 
odic media. Of particular current experimental interest 
are vortex lattices and 2d Wigner crystals [^5[ . In all 
cases, translational and temporal LRO may be destabi- 
lized both by phonons (phase fluctuations) and by topo- 
logical defects (phase slips). Provided reflection invari- 
ance is broken by an external drive field, we expect linear 
gradient terms (e.g. vd x <j>) in the equations of motion. 
Preliminary investigation of driven lattices suggests that 
such terms play a similar role in that case There 
are also examples of oscillatory states in pattern forming 
systems without broken reflection symmetry - leading, 
we expect, to the absence of the ordered phase in the 
presence of disorder and noise. 
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